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1 Introduction 



1.1 Motivation 

Mercury is the inner most of the four terrestrial planets in the Solar system, 
moving with high velocity in the Sun's gravitational field. Only comets and 
asteroids approach the Sun closer at perihelion. This is why Mercury offers 
unique possibilities for testing General Relativity [1] and exploring the limits of 
alternative theories of gravitation with an interesting accuracy [6] . 

As seen from the Earth the precession of Mercury's orbit is measured to 
be 5599.7 seconds of arc per century [4]. As observed, in 1859 by Urbain Jean- 
Joseph Leverrier there is a deviation of Mercury's orbit from Newtonian's predic- 
tions that could not be due to the presence of other planets. Urbain Jean- Joseph 
Leverrier had also been the first to attempt to explain this effect. 

Perihelion precessions of Mercury and other bodies have been the subject of 

experimental study from AD 1765 up to the present. In 1882 Simon Newcomb 

obtained the value 43 seconds per century for the discrepancy for Mercury [5]. 

According to Pireaux et al [6] , the observed advance of the perihelion of Mercury 

that is unexplained by Newtonian planetary perturbations or solar oblatness is 
l 

A u)obs = 42.980 ± 0.002 arc — seconds per century 
2tt(3.31636± 0.00015) x 10" 5 radians 
415.2019rcvolutions 
= 2tt(7.98734 ± 0.00037) x 10~ 8 radians/revolution (1) 

It is the purpose of this paper to provide the first precise calculation of the 
orbit and the perihelion precession of Mercury in General Relativity including 
cosmological constant contributions [2, 11]. We also compare the result with 
observations (1). We assume that the motion of Mercury is a time-like geodesic 
in a Schwarzschild space-time surrounding the Sun. 

Bounds on the magnitude of the cosmological constant from the observed 
value of Mercury's perihelion shift have been estimated using different methods 
of approximation by Islam [12] and Cardona et al [13] with conflicting results. 
The upper bounds obtained were respectively, |A| < 10~ 42 cm~ 2 [12] and |A| < 
10~ 55 cm~ 2 [13]. The latter value compares favourably to the values for the 
magnitute of the cosmological constant obtained from large-scale observations 
[39] and cosmological considerations [11]. 

In [11] a direct connection between general exact solutions of general rela- 
tivity with a cosmological constant in large-scale cosmology and the theory of 
modular forms and elliptic curves [15, 14] was established. 

As we shall see in the main body of the paper, the two-body central orbit 
problem under consideration is very interesting mathematically and becomes 

1 The observations seem to exclude Brans-Dicke theory with ui ~ 5 whose post-Newtonian 
contribution to the perihelion shift would thus have been 39 arc-seconds per century 1) [4, 6]. 
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more involved with the inclusion of the cosmological constant. In particular, the 
integration of the geodesic orbital equations involves the inversion problem of 
Abelian hyperelliptic integrals and the resulting solution which was first studied 
by outstanding mathematicians such as Jacobi, Abel, Weierstrafi, Riemann, 
Gopel, Rosenhain and Baker. The exact treatment of the orbital problem and 
the techniques developed in this paper are general and should be of interest for 
a variety of problems in various fields of cosmology which are discussed in the 
main text. 

The material of this paper is organized as follows. In the rest of the intro- 
duction we review planned satellite experiments that are relevant to the study 
of Mercury and the minor planets. In section 2, starting with the Schwarzschild 
metric, and including the cosmological constant, we derive the time-like geodesic 
equations. In section 3, we solve exactly the geodesic differential equations and 
calculate the precession of the perihelion of Mercury as well as the perihelion 
and aphelion distances of the planet, with and without the cosmological con- 
stant. Our approach involved the solution of Jacobi's inversion problem i) for 
an elliptic integral (case without cosmological constant) and ii) for a genus-2 
hyperelliptic integral (in the general case with A =^ 0) . Here we acknowledge 
the contribution by Whittaker of the exact solution of time-like geodesies in the 
Schwarzschild field with vanishing cosmological constant [40]. However, in his 
account no attempt was made in calculating the physical characteristics of the 
orbit of planet Mercury. For an alternative method of computation of perihelion 
precession in the case of vanishing cosmological constant we refer the reader to 
ref.[16]. In order to make the presentation self-contained we include in the same 
section a comprehensive presentation of the inversion problem. In section four 
we present a discussion and a summary of our main results as well as further 
possible applications. Finally, in an appendix we present the definitions of the 
genus-2 theta functions that solve Jacobi's Inversion problem. 

1.2 Planned Satellite Experiments 

Solar and planetary astronomers and cosmologists are still fascinated with the 
precession of the Perihelion of Mercury and the dynamics and properties of near 
earth objects such as the minor planets (i.e. Icarus). In fact, there is a growing 
interest in this exciting and vibrant field with several satellite missions planned 
for the next few years. Listed below are some of the key satellite experiments 
(see Pireaux et. al. [6] for an overview): 

• GAIA (ESA) 

planned launch 2009, key goal to provide observational data to tackle an 
enormous range of important problems related to to the origin, structure 
and evolutionary history of our Galaxy. The data will create a precise 
three - dimensional map of more than one billion stars throughout our 
Galaxy - a stellar census. GAIA will observe and discover several hun- 
dred thousand minor planets, determine their physical properties and or- 
bital characteristics. Also three earth crossing asteroids, Icarus, Talos 
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and Phaeton, with favourable combination of distance and eccentricity, 
will provide data to more accurately determine the PPN parameters giv- 
ing upper limits on the relativistic and solar quadrupolc contributions to 
the perihelion precessions of solar objects enabling a more stringent test of 
general relativity. The GAIA mission is a hugely ambitious and exciting 
project with a multitude of objectives (see GAIA Study Report [17] for 
details). 

• BepiColombo (ESA) 

planned launch 2010, Mercury Orbitcr Experiment, key objective to study 
the Mercury environment in order to understand the formation and evolu- 
tion of of planets within the solar system, especially the four closest to the 
Sun - Mercury, Venus, Earth and Mars. BepiColombo will measure Mer- 
cury's motion more accurately than ever before, providing data that will 
allow the determination of the precession of the perihelion, thus providing 
one of the most rigorous tests of general relativity [18]. 

• Messenger (NASA) 

planned launch 2004, a Mercury Orbiter experiment whose main goal is 
to study the chacteristics and environment of Mercury in order to better 
understand the fundamental forces that have formed it and the other 
terrestrial planets [19] 

• Picard (C.N.E.S.) 

planned launch end of 2005, main object is to measure the Sun's quadrupolc 
moment contribution to the precession of the Perihelion of Mercury, thereby 
offering another exacting test of general relativity [20] 

Also there are several near Earth Satellites presently operational, LAGEOS, 
LAGEOS II, CHAMP and GRACE. The combined data from these satellites 
will enable the accurate measurement of the perigee advance in the Earth's 
gravitational field (the determination of the PPN parameters 7, /?, in the Earth's 
gravitational field) , providing yet another test of general relativity [7] . 

All of these satellite experiments are expected, to various degrees, to pro- 
vide data that will provide stringent tests for general relativity and to more 
accurately determine solar and planetary geophysical characteristics. 

Our precise results that follow should be of interest to all of the above ex- 
periments. 

2 Schwarzschild solution and time-like geodesies 

Einstein's equations with the cosmological constant A are as follows 

- \g^R = {8ttG/c 4 )T^ + Ag^ (2) 
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where R^ Vl R denotes the Ricci tensor and scalar respectively. Also G denotes 
Newton's gravitational constant and c the velocity of light. For vanishing stress- 
energy momentum tensor T^ v the field equations reduced to 

R^u = A.g M „ (3) 

The motion of a planet according to General Relativity is a time-like geodesic 
in a Schwarzschild space-time [26] surrounding the Sun. The Schwarzschild 
solution for the metric taking into account the cosmological constant (see e.g., 
ref.[42, 12]) is 

ds 2 = c 2 1 5-^ + -Ar 2 eft 2 

c A r 3 

,., 2GM Q 1 2 x 2 

- (1 2-^ + -Ar 2 ) 1 dr 2 

c z r 6 

- r 2 (d9 2 +sm 2 6d<t> 2 ) (4) 

where Mq denotes the mass of the Sun. For zero cosmological constant Eq.(4) 
reduces to the original Schwarzschild metric element [26]. We note at this 
point, that by taking the time-like geodesic equation as a starting point for 
the subsequent calculus represents an approximation to the real situation: it is 
strictly true for point-like bodies with negligible mass. However, since tum <C 
Mq where tum denotes the mass of Mercury the approximation is adequate for 
all practical purposes. 

The geodesic equation has the general form: 

#x* dx«dx? 

It is easily shown that the motion is confined to the equatorial plane = § , 
so that (5) with fj, = 2 is trivially satisfied [42]. Furthermore, we can ignore (5) 
with fj, = 1 in favour of (4), which is a first integral of the geodesic equations. 
The resulting equations can be written as 

dx» dx v 
ds ds 



c 2 



i-^ + 1a^ (!*) 



3 J \ds J 



2 



-r 2 



ds 



(6) 



ds 2 ^ r ds ds ^ ^ 
d 2 t 2(^ffi + lAr-3) drdt _ 

d^ + r(r _o^e i • ----0 (8) 



r{r - 2^P- + ±Ar 3 ) da ds 
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Next, define v := dt/ds and w :— d(f>/ds. Then eqs(7),(8) can be written as 
dw 2 

— + -to = 9 

dr r 



dv 
dr 



1 1 

-- + 



r - + Ur3 



3 

which can be integrated to yield [42] 



(10) 



w — dcfr/ds = (11) 



r 2 

dt 



ds ! _ i^fo + i Ar 2 



(12) 



where L,£ are arbitrary constants. Substituting in (6) we get the following 
equation after definining a new variable u := r _1 

du\ 2 2GM Q 3 2GM Q 1/1 \ /£' 2 -l\ 



c 2 c 2 L 2 3 \u 2 L 2 J \ L 2 

3 Exact solution of the time-like geodesic and 
the perihelion precession 

3.1 Precise calculation of the perihelion advance assuming 
zero Cosmological Constant. 

In the previous section we derived the geodesic equation for the orbital motion 
of the planet Mercury around the Schwarzschild gravitational field of the Sun 
given by eq.(13). In this section we will solve Eq.(13) with a zero Cosmological 
Constant, the non-zero, more general case, will be considered in section 3.3. 
Making the following definitions, 

2GM f) 2GM l-E' 2 
as := —^- t (3 := -^,7 == JJ~ (14) 

our geodesic equation reduces to: 

2 

= a s u° - u +j2-u + j (15) 

Differentiating Eq.(15) with respect to (p and comparing with the Newtonian 
term we get cL = L*m, where Lm is the angular momentum per unit mass of 
the planet [12]. Equation (15) is a cubic equation which can be reduced to the 
Weierstrafi form using the substitution: 

u = —U+^- (16) 
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Then the WeierstraB representation of the geodesic equation becomes 2 




where the WeierstraB cubic invariants are given by: 
1 (2GM Q /c 2 ) 2 2 

= l-i£c 2 
12 4 

1 (2GM Q /c 2 ) 2 2 (2GM /c 2 ) 2 ( E 2 - c 



93 216 48L 2 M C 16 V L 2 M 

= 7^7 - i^ 2 - ^-C(E 2 - c 2 ) (18) 
216 48 16 v ; v 7 

2 

where L := ^-§- and _E = cE'. Then the solution of Eq.(17) is given by: 

U=p(<j> + e) (19) 

and e is a constant of integration. The solution in terms of the original variables 
is given by: 

1 4 , N 1 
«=- = — p(0 + e) + — 20 

r as das 

or 

r= 2£ r 2 i) 

4p(^ + e) + ± ^ ' 

For the calculation of the perihelion precesion and the orbital characteristics 
of Mercury we use the following values for the physical constants: 

c = 299 792 458 m s" 1 , a s = 2.953 250 08 Km (22) 

The data is taken from [8] and [9] , respectively. The value of the speed of light 
in vacuum, given in Eq.(22), is exact, since the meter is the length of the path 
travelled by light in vacuum during a time interval of 1/299 792 458 of a second 
[8]. As free parameters we may use C and E. Then L 2 M — c 2 L 2 = j?a s 3 . Our 
free parameters are mixed through the WeierstraB invariants 32 , 53 with exact 
numbers. 

We will show in the next section that there is a small parameter space for C 
and E that reproduces the observed values for the characteristics of Mercury's 

2 The differential equation that the WeierstraB function satisfies is the equation of an el- 
liptic curve and is given by : (p' (z)) 2 = 4p(z) 3 — g2p(z) — 33, p'(z) = df ^^ , see [11] for 

further details. The inversion of the elliptic integral f U , du = <j>, by the WeierstraB 

\Z4U 3 -g 2 U-g 3 

function, i.e. U = p((j> + e), is a simple case of the problem of inversion for elliptic integrals. 
3 The integration constant L has dimensions of length and thus scales with ag. 
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orbit. Non-physical solutions that predict complex, negative or too small or- 
bital radii and precession values are disregarded. The sign of the discriminant 
A = g\ — Tlg\ determines the roots e*, i — 1, 2, 3 of the elliptic curve: A > 0, 
corresponds to three real roots ej while for A < two roots are complex conju- 
gates and the third is real. In the degenerate case A = 0, (where at least two 
roots coincide) the elliptic curve becomes singular and the solution is not given 
by modular functions. 

We find that the physically acceptable solutions that reproduce the orbital 
data of Mercury correspond to the case where A > 0, e 3 — > — ^,e 2 — > — ^ 
and ei — > 1/6. The analytic expressions for the roots of the elliptic curve 
may be obtained by using the algorithm of Tartaglia and Cardano [43] . Their 
expressions are given in appendix B. The WeierstraB function, p(z), is an even 
meromorphic elliptic function of periods 2u, 2a/ (i.e., p(z + 2uj) = p(z) = p{z + 
2a/), for all complex numbers z). The two half-periods u> and a/ are given by 
the following Abelian integrals (for A > 0) [40]: 



r°° dt , _ . n dt 

J ei V 4<3 - .92* - 93 ' J-oo v/-4t 3 + g 2 t + .93 



(23) 



The value of the WeierstraB function at the half-periods are the three roots of 
the cubic. For positive discriminant A one half-period is real while the second 
is imaginary. The exact expression for the precesion of the perihelion of planet 
Mercury, is given by 

A° TR = 2(o, - tt) (24) 

which is proportional to the deviation of the real half-period uj of the WeierstraB 
modular form, from the transcendental number ir. 

The exact expressions for the minimum distance of planet Mercury from the 
Sun (Perihelion) and its maximum distance (Aphelion) are given by 4 : 

TP = r Perihelion — ~ i~ (25) 

TA = r Aphelion = ~ ~ i~ (26) 

4e 3 + i 

For such orbits the radius vector oscillates between ta and rp, as the argu- 
ment of the WeierstraB function travels along the straight line from ui' to a/ + u> 
and then to a/ + 2a; in the fundamental period region [21]. 

We also note the following. Although by construction the roots of the cu- 
bic are calculated with arbitrary precision, our output for the perihelion and 
aphelion distances can only be displayed with nine significant figures given the 
nine digit accuracy of the Schwarzschild length as- However, their ratio which 
is given by: 

fP = 4e 3 + l/3 

r A 4e 2 + l/3 V 1 



4 We organize the roots as: e\ > e*2 > e-j, 
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constitutes a genuine and precise prediction 5 . 

For a given choice of values for the free parameters, equations (24)-(27) are 
the output of the precise theory assuming zero cosmological constant, for the 
corresponding physical quantities, that should be tested against observations. 



3.2 Theoretical results 

The exact solution of the geodesic equations of the General Theory of Relativity 
outlined in the previous section predicts a set of coupled results, 

Aoj gtr - Perihelion Precession 

rp - distance of closest approach 

rA - distance of farthest approach 

based solely upon the values chosen for the C and E. As we mentioned not 
all of input values for the free parameters provide us with physical solutions, or 
produce realistic data for the orbit of Mercury. In appendix B, we provide some 
examples of solutions that correspond to specific choices of the free parameters 
that do not describe the orbit of Mercury. More specifically, we describe so- 
lutions for the motion of the test particle that are described by special elliptic 
curves for which one of the Weierstrafi invariants, 52, #3, vanishes 6 . In addition, 
we describe choices of the initial conditions that lead to a vanishing discriminant 
A (singular elliptic curves). 

For C = 1.18554647 x lCT 28 cni- 2 s 2 (I? M = 7.35 6 68 x 10 38 cm 4 s" 2 ) and 
E = 0.0299792454176 x 10 12 cm s" 1 we obtain 

e 2 = -0.0833333172749111339, e 3 = -0.083333322753852543388 

. s~"~r d _ 0XCS6C 

A GTR = 4 3 ^ 0017 . 



century 



r P = 4.59766539 x 10 12 cm, 
112 



r A = 6.97872168 x 10 12 cm (28) 

while the two periods are: lu = 3.141592904646, J = 20.40864976 i and the 
period ratio r = 6.496i, and the radius vector oscillates between 4e °;j_i < r < 

j^-t. For C = 1.1848820116975453 x 10~ 28 cm" 2 s 2 , (L 2 M = 7.36080550 x 

4e3 + o 



10 38 cm 4 s~ 2 ) and E = 0.0299792454178 x 10 12 cm s" 1 we obtain, A GTR = 
42.9776-§£fSc. with the half-periods, lu = 3.141592904505435, J = 20.409634 i 

century r ' ' 

and the period ratio t = 6.497i. Also r P = 4.60057668 x 10 12 cm, = 
6.98186816 x 10 12 cm. 

In Tables 1-3 we list the predicted results for A GTR , rp,rA keeping E fixed 
at the indicated value and varying L. We see that by increasing the input value 
for £ one can obtain a slightly larger value for the precession of perihelion. The 

5 In Newtonian theory the orbit of a planet is described by an ellipse with eccentricity e 
and semi- major axis a. For an ellipse, the perihelion (rp) and aphelion (r%) distances, are 

l-r N It N 

a(l — e), a(l + e), respectively, and then e = ^ £ . 

L + r p /r A 

"These elliptic curves have the property of complex multiplication, see [11] for the defini- 
tion. 
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parameters 

C = 

1.1849947026647969 x 10~ 28 cm- 2 s 2 
E = 

0.029979245417779875 x 10 12 cm s" 1 



roots 

0.16666664004116 
-0.083333317282230892 
-0.083333322758930472 



r P 

TA 



predicted results 

A£™ = 42.9817^ 
4.59976206 x 10 12 cm 



6.98207293 x 10 12 cm 



Table 1: Predictions for A^ TH ,rp,rA for the indicated choice for £,E. The 
two half-periods are: uj = 3.14159290452929, J = 20.409059 i and the period 
ratio r = 6.496i. Also L\, = 7.360 1 05 50 x 10 38 cm 4 s~ 2 



parameters 

£ = 

1.1849627128268641 x 10~ 28 cm" 2 s 2 
E = 

0.029979245417779875 x 10 12 cm s" 1 



roots 

0.16666664004188 
-0.08333331728350096 
-0.08333332275837917 



r P 
ta 



predicted results 

A GTR = 4 2.g 8 05M£^£c 

w century 

4.60012605 x 10 12 cm 



6.98170894 x 10 



12, 



Table 2: Predictions for A^ TR ,rp,rA for the indicated choice for £, E. The 
two half-periods arc: u = 3.141592904522524, J = 20.409391 i and the period 
ratio r = 6.4965i. Also L 2 = 7.36 3 4 20 x 10 38 cm 4 s" 2 



parameter E fixes the sum of the radii rp + ta- In fact, in tables 1-3 the sum 
rp + rA is equal to 5.7909175 x 10 12 cm. The latter value is reported as the data 
for the semi-major axis a in ref.[10]. In table 2 the value for Lm that corresponds 
to the choice for C for fixed as is equal to the Newtonian value that would be 
obtained by assuming the orbit is an ellipse with eccentricity e [42, 12]. 

The choice of initial conditions presented in Tables 1-3, lead to a set of 
predictions for the perihelion advance in agreement with observations in Eq.(l). 
In addition, they lead to a set of predictions for the perihelion and aphelion 
distances of the planet, whose half-sum, rp + T ' A , is in agreement with the data 
for the semi-major axis a in ref.[10]. The calculated theoretical values for the 
perihelion and aphelion distances compare favourably with best current values 
for the orbital data for Mercury [27] . More precise measurements of the orbit 
of Mercury as described for instance in the ESA mission BeriColombo, will 
determine the perihelion and aphelion distances as well as the precession of 
the perihelion more accurately and will further restrict the choice of initial 
conditions. 



3.3 Precise calculation of the perihelion advance with the 
contribution of the Cosmological Constant. 

As we saw in section 3.1, the integration of the geodesic equation that describes 
the orbit of Mercury in the central gravitational field of the Sun (assuming 
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parameters 

C = 

1.1848820116975453 x 10~ 28 cm- 2 s 2 
E = 

0.029979245417779875 x 10 12 cm s" 1 



roots 

0.166666640043693405 
-0.083333317286706146 
-0.083333322756987258 



r P 
ta 



predicted results 

A£™ = 42.9776^ 
4.60104489 x 10 12 cm 



6.98079010 x 10 12 cm 



Table 3: Predictions for A^ TR , rp, ta for the indicated choice for £, E. The two 
half-periods are: lu = 3.141592904505435, J = 20.410231661 i and the period 
ratio r = 6.497i. Also L\, = 7.36080550 x 10 38 cm 4 s~ 2 . 



vanishing cosmological constant) involved an elliptic integral, whose inversion 
by the Weierstrafi elliptic function, provided the exact solution in closed analytic 
form for the orbit of the planet eq.(19). 

Elliptic integrals are special cases of the so called Abclian integrals. Accord- 
ing to Jacobi an Abclian integral is an integral of the form J R(x, y)dx, where 
R(x, y) is a rational function in x and y and x, y are connected through an 
equation f(x,y) = 0, where / is an irreducible polynomial. In the special case 
y 2 = P(x), where P(x) is a polynomial of n th degree, with no multiple roots, 
the Abelian integral is called elliptic when n — 3 or n = 4, and hyperclliptic 
when n > 5. 

Including cosmological constant effects we need to calculate the integral: 



du 



2GM Q „3 „2 , 2GM „ 1--E' 2 A A 



or 



f 



u du 



^ 5 _ U 4 + ^ U 3_^ 2 __^_A U2 



(29) 



(30) 



Eq.(29) defines a hyperelliptic integral whose inversion involves genus 2 Abelian- 
Siegelsche modular functions. The problem of inversion (whose solution was 
culminated in the Jacobisches Umkehrtheorem) of hyperelliptic integrals were 
first investigated by Abel [28], Jacobi, Gopel and Rosenhain [29], [30], [31]. The 
explicit solution of Jacobi's inversion problem in terms of higher genus theta 
functions was provided by Gopel and Rosenhain for the case n = 5 or 6, and 
the general solution for the hyperelliptic case (i.e. V n > 5) was provided by 
Weierstrafi [32]. Ricmann introduced the idea of a Ricmann surface to study 
algebraic singularities. He also introduced the Riemann theta function which 
served as a useful tool for solving the Jacobi's inversion problem [33]. For full 
details and extended bibliography to the original literature we refer the reader 
to the book of Baker [34] . 

In what follows we discuss first the Jacobi's inversion problem and then we 
proceed to determine the effect of the cosmological constant on the perihelion 
advance of Mercury as well as its effect on the radii r P , ta- 
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3.4 Abel's theorem and Jacobi's inversion problem 

Let the genus g Riemann hyperelliptic surface be described by the equation: 

y 2 = 4(x - ai) • • • (x - a g )(x - c)(x - d) ■ ■ ■ (x - c g ) (31) 

For g — 2 the above hyperelliptic Riemann algebraic equation reduces to: 

y 2 = A(x — ai){x — a 2 )(x — c){x — C\){x — c 2 ) (32) 

where a\,a 2 , cc\, c 2 denote the finite branch points of the surface. 

The Jacobi's inversion problem involves finding the solutions, for Xi in terms 
of m, for the following system of equations of Abelian integrals [34]: 

U x + ■ ■ ■ + Mi = Ui 

: + ••■ + : : 
u x ^ + ■■■ + u x ^=u g (33) 

where «?•" = f x = f . . . = f x sL^st*. 

1 Ju y ' £ J ft y ' '3 J fj, y 

For 5 = 2 the above system of equations takes the form: 

f Xl dx f X2 dx 

/ — + / — = Hi 

f Xl x dx ^ f X2 x dx _ 

J a\ y J a,2 y 

where u\, u 2 are arbitrary. The solution 7 is given by the five equations [34] 
6 2 {u\u b ' a ) 



e 2 ( u ) 



A(b - xi){b - x 2 ) ■ ■ ■ (b ■ 
A(b - x x ){b - x 2 ) 



= (35 ) 

where /(x) = (x — ai)(x — a2)(x — c)(x — ci)(x — C2), and e 7TlPP — ±1 according 
as u b ' a is an odd or even half-period. Also 6 denotes a finite branch point and 
the branch place a being at infinity [34]. The symbol 6(u\u b ' a ) denotes a genus 2 
theta function with characteristics: 9(u;q,q') [34], where u, = (ui,u 2 ), denotes 
two independent variables, see appendix A for further details. From any 2 of 
these equations, eq.(35), the upper integration bounds x\,x 2 of the system of 
differential equations eq.(34) can be expressed as single valued functions of the 
arbitrary arguments Wi,U2- For instance, 

1 e 2 (u\u a ^ a ) , s 

xi=a 1 + — M— - (36) 

A 1 (x 2 -a 1 ) 2 u) 



7 The definitions of the genus two theta functions appearing in the solution of Jacobi's 
inversion problem are given in the appendix A. 
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and 



x 2 



(0 2 - 0i)(02 + Oi) + ^ V(n) " M^IFW 



2(oi - 2 ) 



± 



where 



/ \/ \ . l 6 2 (u\u a i- a ) 1 e 2 (ti|«°2.») 1 . / n 

{a 2 -a 1 ){a 2 +a 1 ) + — i \l {u) - \\ {u) -4(01-02)77 



?7 := a 2 ai(ai — a 2 ) 



2(ai - a 2 ) 



a 2 6» 2 (u|w ai > a ) ai 6> 2 (u|u a2 '°) 
# 2 ( u ) + A2" # 2 ( u ) 



(37) 
(38) 



Also, Ai = ±- 



V e * ,pp '/'(«.)' 

The solution can be reexpressed in terms of generalized Weierstrafi functions: 



(12) P2, 2 (u) ± ^/p| 2 (u) +4p 2 ,l(w) 
C fc = 9 ' k=l, 



where 



P2, 2 ( U ) = 



/ \ / 1 \ 1 e 2 (u\u a ^- a ) . 1 e 2 (u\u a z- a ) 

(Ol - 2 )(0 2 + Oi) - ^ V(n) + ^ 



01 — a 2 



and 



P2,l( u ) = 



-Oi0 2 (oi - o 2 ) - ^ &A{u) + — x 9Hu) 



01 — a 2 



(39) 



(40) 



(41) 



Thus, xi,x 2 , that solve Jacobi's inversion problem Eq.(34), are solutions of a 
quadratic equation [29, 34] 



Ux 2 - U'x + U" = 



(42) 



where U, U', U" are functions of ui,u 2 . In the particular case that the coefficient 
of x 5 in the quintic polynomial is equal to 4, U = 1, U' = p 2 , 2 (u), U" — p 2 ^(u). 



3.5 Inversion of the hyperelliptic integral. 

Let us apply the solution of the Jacobi's inversion problem given by eqs.(35)- 
(37) to the hyperelliptic integral (30). We have the following correspondence 
between the variables of the hyperelliptic integral and the variables in (34): 
the upper integration bounds in (34) Xi, correspond to Ui and the periods Ui 
correspond to <p im Then for: 

0i = 4>r 2 

4> 2 = r 2 a2 (43) 
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u\ = ai, U2 — u, u arbitrary. (44) 

Thus we have that: 



1 6» 2 ($|$ 01 >°) 
~M 6> 2 ($) 

1 6» 2 ($|$ a2 ' a ) 
~M 2 ($) 



= (ui - ai)(u 2 - ai) = 

= (en - a 2 )(u - a 2 ) (45) 



and 



and the solution for u in the hyperelliptic integral (j> 2 — </> 2 ' a2 — J" 2 (where 
y is the radical of the quintic polynomial 8 ) is given by 

«2 = u = — - — + (46) 

where 

1 6» 2 (<J>|$ a2 ' a ) 

P2^)=ai+a 2 + — r v ' - (47) 

A 2 (oi-a 2 ) 

«i <9 2 ($|$ a2 > a ) 

p 2 i W = -aia 2 — - — -. — -— ^ '- (48) 

V ; A 2 (ai-a 2 ) 2 ($) v ' 

3.6 Roots of the quintic, periods of the hyperelliptic sur- 
face and the cosmological constant effect 

In principle, the sign of the cosmological constant is an additional parameter, 
besides its magnitude. Recent observations of large-scale cosmology indicate a 
positive cosmological constant of magnitude <~ lCT 56 cm~ 2 [39]. 

In the presence of the cosmological term there are five branch points for 
the hyperelliptic surface, eq.(32), which are obtained by solving the quintic 
polynomial that appears in the time-like geodesies, Eq.( 29). 

For negative cosmological constant all the roots are real. For positive cosmo- 
logical constant and depending on its magnitude and the values of the parame- 
ters £, E three roots are real and two complex conjugates. For some particular 
values all the roots are real. 

Let us start our discussion with the case of the negative sign. Let us arrange 
the roots in ascending order of magnitude and denote them by e2 9 , e 2g _i, • • • , eo, 
2, so that e 2 j,e 2 j_i, are respectively c s _j+i, a s _j+i and eo is c. We then de- 
fine linearly independent Abelian integrals of the first kind [34], denoted by 
Z7f' a , i = l,---g, whose periods we want to calculate. These integrals are 
such that dll^' a /dx — tp r /y, where ip r is an integral polynomial in x, of degree 
g — 1 — 1 at most, with only real coefficients. Then the half-periods, l^ 4 ' 63 and 
Z7® 2 ' ei , are real, while the half-periods C/^ 3,e2 and U^ 1,e ° are purely imaginary 

8 we can always bring the quintic polynomial of the hyperelliptic integral into the defining 
form (32) of the genus-2 hyperelliptic surface. 
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[34]. For clarity, by we denote Q U^ es denotes Q ^ and similarly 

for the rest of the periods. 

Every other period can be expressed as a linear combination of the above 
Abelian integrals. For instance, the period U^ 1,a is given by the equation: 

U^' a = -U^ 2 ' ei - C/ r e4 < es + C/ r ei ' eo (49) 

while the period U^°' a — U^"'°° is given by 

Jje M = ^64 ,e 3 ( 50 ) 

For A = -10- 55 cm- 2 ,£ = 1.1848820116975453 x 1(T 28 cm" 2 s 2 (L 2 M = 
7.36 8 5 5 x 10 38 cm 4 s" 2 ), E = 0.0299792454178 x 10 12 cm s" 1 the roots are: 

e 4 = -1.143376827586711663025 x 10" 24 

e 3 = 1.14337682758671162302505 x 10" 24 

e 2 = 1.43228141354887622443346 x 10" 13 

d = 2.173640544970253310548427 x 10" 13 

e = 0.000003386099606939 (51) 



The half-periods are calculated to be: 

r e 3 ,e 2 

C/ 2 ei ' e2 = 3.14159290450544 (52) 



C/ 2 e3 ' e2 = 2.264048243121482 i, U^ ue ° = 20.4094269 i 



the real half period: f 61 ^ = 3.14159290450544, leads to Aff™ = 42.9776 arc 7 scc 

L J e2 y ' ^ century 

in agreement with observations (1). 

For A = -10- 42 cm- 2 ,£ = 1.1848820116975453 x 10" 28 cm" 2 s 2 ,E = 
0.0299792454178 x 10 12 cm s" 1 the roots are: 



e 4 


= -3.615599293152843540392 x 10 


-18 


e3 


= 3.6157507116827139496388 x 10" 


18 


e2 


= 1.4322814108727319240989 x 10" 


13 


ei 


= 2.1736405461322123121789 x 10" 


13 


eo 


= 0.000003386099606939 





(53) 



The periods are calculated to be: 

J7 2 e3 ' e2 = 2.26404824460 i, J7 2 ei ' e ° = 20.4094269 i 

U* ue2 = 3.14159290255, C/ 2 e3,e4 e R (54) 

Here there is substantial effect on the perihelion advance due to the cosmological 
constant with A G , TR = 42.6427 arc 7 scc in conflict with observations (1). 

& century v 7 
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A = 


A = -10- 42 cm- 2 


A = 10- 42 cm- 2 




3.141592904505435 


3.14159290255 


3.141592903870775 


\GTH 

'-^LlJ 


42.9776 ar ? oc 

century 


42.6427 arc f soc 

century 


42.8689 arc f soc 

century 


r P 


4.60057668 x 10 12 cm 


4.60057668 x 10 i2 cm 


4.60057668 x 10 i2 cm 


T A 


6.98186816 x 10 12 cm 


6.98186817 x 10 12 cm 


6.98186815 x 10 12 cm 



Table 4: Theoretical predictions for perihelion precession for A ^ 0, and C = 
1.1848820116975453 x 10~ 28 cm- 2 s 2 , E = 0.0299792454178 x 10 12 cm s" 1 . For 
comparison we list the A = case for the same set of values for C, E. For this 
choice of C, L 2 M = 7.36080550 x 10 38 cm 4 S - 2 . 



For A = 10- 42 cm- 2 , C = 1.1848820116975453x 10~ 28 cm" 2 s 2 ,E = 0.0299792454178x 
10 12 cm s _1 the roots are: 

ei = 0.000000010574713017622 x — 

a s 

e 2 = 0.0000000160482602747 x — 

e 3 = 0.2499999733770 x — 

a s 

e 4 = (-5.58970979303625 x 10~ 18 - 2.6694981178048 x 10~ 13 i) x — 

a s 

e 5 = e 4 (55) 
with the periods f ei ^ = 3.141592903870775, f 62 ^ = 20.4095744i etc. In 

Je2 y ' J e% y 

this case we obtain 2tt(7.96669 x 10" 8 )radians/revolution = 42.8689 arc 7 scc in 

v ' 1 century 

conflict with observations. For A = 10- 56 cnr 2 ,£ = 1.1848820116975453 x 
10~ 28 cm- 2 s 2 cm 4 s~ 2 , E = 0.0299792454178 x 10 12 cm s" 1 the roots arc: 

ei = 0.0000000105747129978 x — 

e 2 = 0.0000000160482602833 x — 

a s 

e 3 = 0.2499999733770 x — 

a s 

e 4 = (-5.5897098056 x 10~ 32 - 2.669498120170 x 10~ 20 i) x — 

e 5 = e 4 (56) 

with the real half period: Q ^ = 3.141592904505435, which leads to A^ TR = 
42.9776 arc 7 scc in agreement with observations (1). 

century v / 

We summarise our results in Table 4 and 5. We also repeated our calculations 
for the values for C,E as in Tables 1 and 2, and we found: For the values 
for C,E as in table 1 for A = -lO^cm" 2 , f ei ^ = 3.1415929045292983, 
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A = -10- 55 cm- 2 


A = 10- 5B cm- 2 


A = 




3.14159290450544 


3.141592904505435 


3.141592904505435 


\GTH 

'-^LlJ 


42.9776 ar ? oc 

century 


42.9776 ai T° c 

century 


42.9776 arc f soc 

century 


r P 


4.60057668 x 10 12 cm 


4.60057668 x 10 i2 cm 


4.60057668 x 10 i2 cm 




6.98186816 x 10 12 cm 


6.98186816 x 10 12 cm 


6.98186816 x 10 12 cm 



Table 5: Theoretical predictions for perihelion precession for A ^ 0, and C = 
1.1848820116975453 x 10~ 28 cm- 2 s 2 , E = 0.0299792454178 x 10 12 cm s" 1 . For 
comparison we list the corresponding results for the A = case. For this choice 
of C, L 2 M = 7.36080550 x 10 38 cm 4 S - 2 . 



a gtr. = 42. 98 i7^£_^ an d for A = 10- 56 cm- 2 , f ei ^ = 3.141592904529299 

<*> century ' Je2 y 

and A2 TR = 42.9817 arc 7 scc . 

w century 

For the values for C,E chosen in table 2, for A = -lO^cm" 2 , f 1 ^ = 

' J €2 y 

3.14159290452253, A% TR = 42.9805 a ^~ s r e y c , while for A = 10- 56 cnr 2 , Q ^ = 
3.141592904524534 and A% TR = 42.9809 ^~ s r ° y c . Thus we see that the cosmo- 
logical constant depending on its magnitude has an effect on the perihelion 
precession and in fact is in conflict with observations for magnitudes |A| ~ 
10 2 cm -2 and larger while current favoured values of the cosmological con- 
stant [39, 36, 37] are compatible with observations of the Mercury's perihelion 
advance and do not have significant effect either on the perihelion precession or 
the radii rp, ta- Evidently, more precise observations can put a rigorous upper 
bound on the magnitude of the cosmological constant. 



4 Discussion and further applications 

In this paper we dealt with the motion of a test particle in a central gravita- 
tional field. In particular we investigated the orbit of planet Mercury around 
the gravitational Schwarzschild field of the Sun by solving exactly the geodesic 
differential equations with and without the cosmological constant and deter- 
mined the perihelion advance. As we saw in the main body of the paper, the 
analytic solution of the orbital problem is provided by the solution of the inver- 
sion problem first enunciated by Abel and Jacobi in their study of elliptic and 
hyperclliptic modular functions. 

The perihelion precesion of the planet in the exact solutions obtained de- 
pends on fundamental properties of the modular functions and in particular 
their periods. 

We have also calculated the precise perihelion and aphelion distances of the 
planet. The calculated theoretical values compare favourably with best current 
values for the orbital data for Mercury [27]. 

For zero cosmological constant the exact solution for the orbit of Mercury 
is expressed by the Weierstrafi Jacobi modular form. The orbital data of the 
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planet are reproduced for positive discriminant (A > 0) and when all roots 
of the cubic are real. For the particular choice of the free parameters given in 
Tables 1-3 the theoretical predictions for the perihelion advance are in agreement 
with observations Eq.(l). The complex structure of the torus that describe the 
particular orbital solutions is r = 6.497 i. We note that the modular properties 
of the exact solutions in large-scale cosmology including cosmological constant 
effects allow for both signs of the discriminant A and in some interesting cases 
the complex structure of the corresponding tori is a fixed point of the modular 
group [11]. 

For non-zero cosmological constant the exact solution for the orbit is pro- 
vided by quotients of genus two theta functions and therefore is described by a 
hypcrclliptic curve of genus two. Alternatively, the solution can be reexpressed 
in terms of generalized Weierstrafi functions. We investigated the effect of var- 
ious values for the cosmological constant on the perihelion precession. Mag- 
nitudes of A <~ 10~ 42 cto~ 2 arc in conflict with observations (1), while current 
favoured values from large scale cosmology are compatible with observations. 
More accurate measurements will place a rigorous upper bound on the magni- 
tude of the cosmological constant. 

The results of this paper, are in full agreement with the method suggested in 
[11] namely: the properties of modular theta functions associated with Riemann 
surfaces and in particular the non-linear differential equations they obey, provide 
a new way to deal with Einstein' equations. In this way, the non-linearity of 
the gravitational field equations can be tackled effectively and general exact 
solutions be generated. 

Our exact treatment of the geodesic equations can be applied to a variety 
of problems in cosmology where the two body central orbit problem in general 
relativity is encountered. More specifically our techniques can be applied to 
address problems such as : the determination of orbital parameters of stars 
around massive gravitational objects such as Neutron Stars [22] or candidate 
Black Holes [23], galactic dynamics including the missing mass problem [3] which 
will be a subject of separate publication , or other scales in cosmology [24, 38]. It 
would also be interesting to generalise the techniques developed in this paper to 
the gravitational field of a rotating mass. Thus effects, like the Lense-Thirring 
[25] precession associated with such gravitational fields could be investigated. 
We note that the Schwarzschild solution because of its static character does not 
describe rotation of the mass distribution. Such generalisation is of additional 
interest for the following reason: taking into account very fine corrections due 
to a non-zero cosmological constant calculated in this paper, one should look 
carefully if they are not comparable with other subtle relativistic effects, like 
the rotation of the Sun and the associated Lense-Thirring effect. Other subtle 
effects include the constant loss of mass during cosmological times. 

In addition, the precise general relativistic predictions obtained can be com- 
pared with the predictions of scalar-tensor type gravity theories such as the 
Brans-Dicke theory. We already mentioned that the latter theory for small val- 
ues of its additional parameter uj (u <~ 5) is in conflict with the observations 
for the perihelion precession of Mercury Eq.(l) and current estimations of the 
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Sun's quadrupolc moment [6, 4]. 

Non-commutative field theories are also strongly constrained from the pre- 
cession of the perihelion of Mercury [35]. 

We are entering a new era in cosmology where precise theoretical results can 
be compared with the ever increasing accuracy of observations. 
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A Definitions of genus-2 theta functions that 
solve Jacobi's inversion problem 

Riemann's theta function [33] for genus g is defined as follows: 

Q(u) := e 2 * mn+l7inn2 (57) 

m,- ,n g 

where fin 2 :~ flun\ + ■ ■ -2fli 2 nin2 + • • • and un :— u\n\ + ■ ■ -u g n g . The 
symmetric g x g complex matrix fl whose imaginary part is positive definite is 
a member of the set called Siegel upper-half-space denoted as Cs ■ It is clearly 
the generalization of the ratio of half-periods r in the genus g = 1 case. For 
genus g = 2 the Riemann theta function can be written in matrix form: 



e(u,fi) = e mtnnn+27Titi 



n£Z 2 



E 



7ri ni n2 



flu fli2 

fll2 fl22 



7ll 

n 2 



+2iti( m n 2 



Ul 

u 2 



Riemann's theta function with characteristics is defined by: 

Q(u:q,q') := ^ e 2-Kiu(n+q')+iTTn(n+q') 2 +2lTiq(n+q) 



(58) 



(59) 



herein q denotes the set of g quantities qi , ■ ■ ■ ,q g and q' denotes the set of g 
quantities q[,--- ,q' g . Eq.(59) can be rewritten in a suggestive matrix form: 



e 



q' 

q 



(u,fl)= Y e* it ( n+q ' )n{n+q ' )+2mt{n+q ' ){u+q \ q,q' £ 



(60) 



nezs 
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The theta functions whose quotients provide a solution to Abel-Jacobi's 
inversion problem are defined as follows [34]: 

9(u;q,q') := ^ e a « 2 + 2 M™+9')+b(™+</) 2 +2'™<j(n+9') (51) 

where the summation extends to all positive and negative integer values of 
the g integers m, • • • , n g , a is any symmetrical matrix whatever of g rows and 
columns, h is any matrix whatever of g rows and columns, in general not sym- 
metrical, b is any symmetrical matrix whatever of g rows and columns, such 
that the real part of the quadratic form bm 2 is necessarily negative for all real 
values of the quantities mi, • • • , m g , other than zero, and q,q' constitute the 
characteristics of the function. The matrix b depends on \g{g + 1) independent 
constants; if we put iiril = b and denote the (/-quantities hu by inU, we obtain 
the relation with Ricmann's theta function: 

e(u;q,q')=e au2 e(U;q,q') (62) 

The dependence of genus-2 theta functions on two complex variables is de- 
noted by: 9(u;q,q') = 9{u\,U2',q,q'), the dependence on the Siegel moduli 
matrix Q by: 9{u\,U2, ^; q, q')- To every half-period one can associate a set of 

characteristics. For instance, the period u a - ai = i ^ j jj ^ while 9{u) is a 

theta function of two variables with zero characteristics, i.e. 9{u) — 9(u; 0, 0) = 

q |j {u, Q). Also, Weierstrafi had associated a symbol for each of the six 

odd theta functions with characteristics and the ten even theta functions of 
genus two. For example, 9(u) is associated with the Weierstrafi symbol 5 or 
occasionaly the number appears as a subscript, i.e. 9(u) 5 . 

The matrix elements hij , fi ^ can be explicitly written in terms of the half- 
periods Ur' a defined in section 3.6. For instance, the matrix element hn = 

x in, while Sin 



Let us perform some consistency checks for particular cases of the inversion 
problem: 

du 

— = <Pi 
V 

-^0 2 (63) 
V 

For instance for 4> 1 = 0^ 2 ' ai ,0 2 = <j)£' ai , U2 =02 = ei,ui = C2 = Let 
us check if the formulas eqs.(35)-(37) reproduce the chosen upper limits of the 
period integrals in the inversion problem. 
In this case, 

1 9 2 (^ ai - a ) 



r ui du 


fM2 


J ai y 




1 u du 


r U2 u 


y 


J 0,2 


a i J, 




A2 = 





A x 9 2 {cj ) ) 
1 9 2 {cj>\^< a ) 
A 2 9 2 ((f>) 



(ui - ai)(u 2 - ai) = (c 2 - fli)(a 2 - 01) 
(64) 
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and therefore using (37) we obtain u 2 —0,2- Then from 

1 2 (0|0 ai < Q ) 1 . 

ui = 01 + — — -27— — = ai H (c 2 -oi)(a 2 - ai) = c 2 . (65) 

^i(ti2-ai) fl 2 ^) a 2 -ai 

B Roots of the cubic and special orbit cases for 
vanishing cosmological constant 

The roots of the elliptic curve can be calculated analytically using the algorithm 
developed by Tartaglia and Cardano [43]. Their general expressions are given 
by: 



n 



T2 = 



2 1 / 3 1 2 + 18£c 2 - 27CE 2 + 3V3VJ j 2/3 + 2 - 6£c 



12 2 2 /3 I (2 + 18£c 2 - 27CE 2 ) + 3^/3^} 

r 1 2/3 

p2V3 1 2 + 18£c 2 - 27CE 2 + 3^/3^} + p 2 (2 _ 6£c2) 

1/3 

12 2 2 /3 {(2 + 18£c 2 - 27CE 2 ) + 3^/3^} 

r i 2/3 

p 2 2V3 1 2 + I8£c 2 - 27CE 2 + 3^3^} + p(2 - 6Cc 2 ) 

i - To 

12 22/3 {(2 + 18£c 2 - 27CE 2 ) + 3^/3^} 



where 



S := C (8c 4 £ + 4C 2 c 6 + c 2 (4 - 36E 2 C) + E 2 (27E 2 C - 4)) (67) 
and p = e 2m / 3 . 

Equation (21), represents the exact solution in closed analytic form for the 
orbit of a test particle in the central gravitational field of the Sun assuming zero 
cosmological constant. As it was mentioned in section 3.1, not all choices of 
initial conditions will lead to planetary orbits. 

Below we list the orbits that correspond to specific choices of initial condi- 
tions and in particular to the two elliptic curves for which one of their Weierstrafi 
invariants vanishes. 

i) Case of negative discriminant and g 2 — 0,53 7^ 0. 

In this case for C = or equivalently, L 2 M — 3c 2 a|, we obtain: g 2 — 
0, 93 = 5jg — j^z — j<^z(E 2 — c 2 ). The discriminant A = — 27<?f < and two 
roots are complex conjugates and one is real. The Weierstrafi function p((f>) is 
real along the diagonal of the fundamental period parallelogram (FPP), and the 
real root e 2 is located at the point ui + <J . As the argument of the Weierstrafi 
function changes from ui + w' to along the diagonal of FPP, its value changes 
from e 2 to 00. 

For E = ^/|c, also 53 = and then there is no elliptic curve since the 
discriminant vanishes and all three roots are zero. 
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ii) Case of g% = 0, 52 > 0. 

In this case there are three real roots, one of them e 2 = and the other two 
are equal in magnitude and opposite in sign, e\ = ^ 1 ^^ c , e 3 = —e\. Here, 

as the argument of Weierstrafi function p, travels along the line from u> + u>' 
to lu in the FPP, its value changes from to e\ > 0. Thus, the radius vector 
varies from r e2 = 3as to r ei = a * j . On the other hand in the region, of the 

4ei + 3 

FPP along the line from lj + lj' to lj' + 2lj the Weierstrafi function is negative 
e 3 < p{4>) ^ 0- The radius vector increases from 3ag and tends asymptotically 
to infinity when 4p(</>) + 3 — > 0. 

iii) Case of vanishing discriminant. 

The discriminant A = g\ — 11g\ can be written as follows: 

£(-4£ 2 c 6 + £(-27£ 4 + 36c 2 £ 2 -8c 4 )+4(i; 2 -c 2 )) , 

A = — (68) 

256 v ' 

It vanishes, besides the case we mentioned above for which g 2 = .93 = 0, also 
for £ = and for 

£ = (8c 4 - 36c 2 S 2 + 27£ 4 ± ^-BUc^E 2 + 1728c 4 ^ 4 - 1944c 2 £ 6 + 729^ 8 ) 
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